Molecular Adsorption on Metal Surfaces with a van der Waals Density Functional 



Guo Li, 1 ' 2 ' 3 Isaac Tamblyn, 4 Valentino R. Cooper, 5 Hong-Jun Gao, 2 and Jeffrey B. Neaton 4 

1 International Center for Quantum Design of Functional Materials (ICQD)/Hefei 
National Laboratory for Physical Sciences at the Microscale (HFNL), 
University of Science and Technology of China, Hefei, Anhui, 230026, China 
2 Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
3 Department of Physics and Astronomy, University of Tennesse, Knoxville, TN 37996, USA 
4 Molecular Foundry, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 
5 Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 

The adsorption of 1,4-benzenediamine (BDA) on the Au(lll) surface and azobenzene on the 
Ag(lll) surface is investigated using density functional theory (DFT) with a non-local density func- 
tional (vdW-DF) and a semi-local Perdew-Burke-Ernzerhof (PBE) functional. For BDA on Au(lll), 
the inclusion of London dispersion interactions not only dramatically enhances the molecule- 
substrate binding, resulting in adsorption energies consistent with experimental results, but also 
significantly alters the BDA binding geometry. For azobenzene on Ag(lll), the vdW-DF produces 
superior adsorption energies compared to those obtained with other dispersion corrected DFT ap- 
proaches. These results provide evidence for the applicability of the vdW-DF method and serves as 
a practical benchmark for the investigation of molecules adsorbed on noble metal surfaces. 



Understanding the fundamental interactions that bind 
organic molecules to noble metal substrates is of cru- 
cial importance in molecular-scale electronics and self- 
assembly, where the competition between molecule- 
substrate and intermolecular interactions can lead to 
templated arrangements with specific spectroscopic and 
transport properties.^1 As the forces driving the for- 
mation of these organic-inorganic assemblies often in- 
clude both local chemical bonding and non-specific long- 
range interactions, it is essential to have an accurate de- 
scription of both contributions (e.g. Ref[3). While den- 
sity functional theory (DFT) provides a many-particle 
framework that, in principle, incorporates both local 
and nonlocal interactions, common semi-local approxi- 
mations to DFT neglect long-range attractive contribu- 
tions to van der Waals interactions, so-called "London 
dispersion forces". 45 However, in recent years, progress 
has been made towards including London dispersion cor- 
rections within standard DFT. These approaches run the 
gamut from semi-empirical methods to the development 
of more accurate exchange-correlation functionalsP^ 
Among these methods, a fully first-principles van der 
Waals density functional (vdW-DF)PH2] has been devel- 
oped to accurately include the effects of London disper- 
sion forces. This method has been shown to be rela- 
tively accurate as well as computationally tractable and, 
as such, has been applied with success to a wide range of 
systems including the adsorption and wetting of various 
sur faces! EHHI 

In the present work, we perform DFT calcula- 
tions, with and without long-range London disper- 
sion corrections, to investigate the adsorption of 1,4- 
benzenediamine (BDA) and azobenzene on the (111) sur- 
faces of Au and Ag. In the BDA/Au(lll) system, we 
demonstrate that inclusion of dispersion forces via the 
recent vdW-DF results in enhancements to molecule- 
substrate binding, bringing predictions for adsorption en- 
ergies into agreement with experimental results.^ Fur- 



thermore, use of the vdW-DF significantly alters the pre- 
ferred orientation of the molecule relative to the Au(lll) 
surface, resulting in a preference for a flat adsorption ge- 
ometry over the tilted configuration previously obtained 
with a generalized gradient approximation (GGA)J^ 
Likewise, in the azobenzene/ Ag(l 11) system, vdW-DF 
results in better adsorption energies than those obtained 
with other semi-empirical dispersion corrections. Com- 
pared with experiment 16 and complementing previous 
vdW-DF studies of flat molecules on Au(lll) surfacespl 
our results provide evidence for the utility of vdW-DF 
for studies of molecule-metal binding. 

Our DFT calculations used both the vdW-DF^ and 
the Perdew-Burke-Ernzerhof (PBE]P GGA to compare 
the effects of dispersion interactions on the adsorption 
of molecules to noble metal surfaces. All calculations 
of BDA molecules on Au(lll) were performed using a 
408 eV planewave cutoff and ultrasoft pseudopotentialJ^ 
as implemented in a modified version of the Quantum 
Espresso simulation package (QE ver. 4. 2.1). 20 (The 
vdW-DF module was obtained from the SIESTA simu- 
lation package. 10 ) Calculations of azobenzene molecules 
on Ag(lll) were performed with the VASP (5.2.12) sim- 
ulation package, employing a 500 eV plane- wave cut off 
and PAW potentials! 21 * 2 ^ For both systems, a periodic 
4-atom-layer slab with 20 A of vacuum and a 2x2x1 
Monkhorst-Pack Appoint mesh were used; both are found 
to result in converged energetics and binding geometries. 
During relaxations, the bottom two layers are fixed and 
all other atoms were allowed to relax unconstrained until 
the forces on each atom were less than 3 me v/A. 

Our computed Au lattice constants are 4.14 A (PBE), 
and 4.25 A (vdW-DF), respectively. The overestimate 
of the lattice constant relative to experiment (4.08 A) 
within vdW-DF has been noted before^HUl anc [ was at- 
tributed to excessive exchange (resulting in unphysically 
strong short range repulsion). We find similar agreement 
with experiment for the bulk Ag lattice constant yielding 



values of 4.16 A, and 4.26 A within the PBE and vdW- 
DF. To model the Au substrate, a 4-layer 4x4 in-plane 
unit cell containing 64 Au atoms was used. Following 
Ref. [I6j we used a 4-layer Ag slab with a 3x6 in-plane 
unit cell of 72 Ag atoms. 

BDA-Au junctions have been intensely investigated 
as a prototype for understanding charge transport at 
the molecular scale ll^MMl i n a BDA-Au junction, 
amine groups preferentially bind to under-coordinated 
Au atoms, resulting in well-defined conduct ance JMUMQJ 
Temperature-dependent helium atom scattering experi- 
ments reported the binding energy (BE) of a related but 
distinct system, sparse BDA sub-monolayers adsorbed on 
flat Au(lll) substrates, as roughly 1 eVp^ stronger than 
the computed bond strength between an amin e group 
and an under-coordinated Au atom (0.4-0.7 eV)! 26 | 2T | ^ 
Moreover, flat Au(lll) is expected to be chemically inert, 
and thus a primary contributor to the large adsorption 
energy is expected to be London dispersion interactions 
between the BDA molecule and the Au substrate. These 
interactions were not explicitly accounted for in previous 
calculations (see Ref. 15). In what follows, we use cal- 
culations with the vdW-DF and the PBE functional to 
compare their performance with experiment, and to bet- 
ter understand the role of nonlocal, dispersion forces on 
BDA adsorption. 

The amine groups at either end of the gas-phase BDA 
molecule adopt a pyramidal structure with two H atoms 
located on one side of the phenyl plane and an electron 
lone pair on the other. Consequently, BDA is stable in 
both the trans- or in the cis- structure (see Fig. Ilk). 
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FIG. 1. (Color online), (a) The trans- and cis- conforma- 
tion of BDA. (b) Dependence of adsorption energy on the 
tilt angle, a, within PBE (black ■) and vdW-DF (red A) for 
trans-BDA on Au(lll). 

Adsorption energies of a trans-BDA on Au(lll) as a 
function of tilt angle relative to the surface are shown in 
Fig. [TJ), where a denotes the angle between the phenyl 
plane of BDA and the Au surface. For each angle, each of 
the atoms within the BDA molecule is constrained to pre- 
serve the tilt angle with the surface. Using this approach, 
the PBE energetic minimum is found to be ~ 35° with a 
binding energy of 0.37 eV. Conversely, the vdW-DF cal- 
culations exhibited significant increases in the adsorption 
energies as well as a shift in the tilt angle to 15°. Ev- 
idently the inclusion of nonlocal dispersion interactions 
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FIG. 2. (Color online). Equilibrium tilting configurations 
obtained within (a) PBE and (b) vdW-DF scheme for the 
trans conformation of BDA on the Au(lll) surface. 



results in the molecule being effectively pulled closer to 
the Au surface. 

To obtain accurate equilibrium adsorption energies and 
optimized configurations, full structural relaxations are 
also performed. For PBE (see Fig.J^a)), we find that the 
adsorption energy is 0.41 eV with a N-Au surface atom 
distance of 2.55 A (at a tilt angle of 29°). These results 
agree with those reported previously, 15 where the binding 
was attributed to a weak but non-negligible amine- Au 
bond. On the other hand, for vdW-DF (see Fig. ^p) the 
adsorption energy is 0.94 eV, in good agreement with 
the experimental value of 1 eV. Here the N-Au distance 
is 3.12 A significantly larger than the PBE result. For 
the relaxed trans molecule on the Au(lll) surface, we 
find that the tilt angle is now 8° with the four amine H 
atoms located at the same height relative to the surface. 
This indicates that the phenyl plane tilt is attributable to 
the £rans-structure rather than the amine- Au interaction 
(see Fig.^). 

BDA can also bind to the Au(lll) surface (see Fig.|3^) 
in the cis-conformation. In this case, we find that PBE 
and vdW-DF give similar lowest energy geometries, with 
the four N-bonded H atoms pointing towards the sur- 
face. The dependence of adsorption energy on the BDA- 
Au(lll) separation is illustrated in Fig.[3]3. PBE calcula- 
tions find an optimum separation distance of 3.71 A with 
a binding energy of 0.39 eV, similar to the tilted struc- 
ture discussed above. vdW-DF calculations predict an 
optimal separation of 3.57 A , and an adsorption energy 
of 0.98 eV, slightly stronger than that of the vdW-DF 
tilted £rans-configuration. The similar binding energies 
between the cis- and trans- conformations obtained us- 
ing PBE, leading to meta-stable tilts, emphasize the fact 
that the molecule-surface interactions are dominated by 
dispersion forces, which stabilize flat adsorption geome- 
tries. Furthermore, the large enhancement in binding 
energies within vdW-DF (0.60 eV) further confirms the 
need for an accurate treatment of non-local interactions 
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FIG. 3. (Color online). Computed binding energy curves 
for the ds-BDA on the Au(lll) surface. The vdW-DF (SC) 
results are computed using a fully self-consistent approach at 
the vdW-DF Au lattice constants. vdW-DF (non-SC) results 
are performed using the PBE charge density (structures and 
lattice constants are based on PBE), i.e. post-processing of 
vdW-DF. 



also the dispersion-corrected methods of Grimme^ and 
Tkatchenko- Schemei^. (These prior results are repro- 
duced in Fig. [4] for comparison). 

Using vdW-DF, we compute an azobenzene-Ag(lll) 
BE of 0.98 eV, in excellent agreement with experiment. 
This is in stark contrast with the small BE predicted 
by PBE (0.1 eV). The larger discrepancy between PBE 
and vdW-DF BE for azobenzene (compared with BDA) is 
consistent with its two polarizable aromatic rings (rather 
than just one). Like PBE however, the vdW-DF pre- 
dicted adsorption height is too large (by about 0.5 A). 
Although this is an improvement over PBE (a 26% over- 
estimate), the adsorption height is still an order of mag- 
nitude larger than typical errors associated with DFT 
bond lengths. This is consistent with vdW-DF's typical 
overestimation of separation distances (a consequence of 
excessively repulsive exchange interaction s 38 * 40 * ^ ) . Our 
computed value for uo (see Fig. [4| is also consistent with 
the experimental measurement. As with BDA on Au, the 
site dependence of the BE seems to be small. Results are 
nearly identical when one of the N atoms is placed above 
an atop site (dN-Ag = 3.63 and BE=0.97 eV). 

Since both the adsorption height and angle have been 
measured experimentally, we can also examine how well 
vdW-DF describes the interactions between the aromatic 
ring and the surface at fixed d^Ag- When we constrain 



in investigations of metal- molecule interfaces. 

Consideration of the adsorption of cis-BDA at dif- 
ferent sites on the Au(lll) surface indicates that the 
Au(lll) surface is energetically fiat for BDA molecules, 
i.e. the adsorption energies for different sites are simi- 
lar (< 10 meV). Although the ordering of binding sites 
within vdW-DF differs from PBE, in both cases the dif- 
ferences between sites is very small, essentially on the 
order of the expected DFT error. Similar results are ex- 
pected for the £rans-configuration. 

To compare with our BDA-Au(lll) calculations and 
further assess the efficacy of the vdW-DF, we also 
compute the equilibrium adsorption geometry and BE 
for a related molecule, azobenzene, on the Ag(lll) 
surface. This system is chosen not only because its 
optical-structural properties have attracted significant 
attention,^ but also due to the recent availability of 
experimental measurements of its adsorption energetics. 
Mercurio et al 16 reported both the binding geometry 
and adsorption energy of azobenzene on Ag(lll) using 
a normal-incidence X-ray standing wave (NIXSW) ap- 
proach. They find that upon adsorption the N-Ag dis- 
tance is 3.07 ± 0.02 A and that the molecule lies fiat 
relative to the surface (uj — — 1±0.2). Furthermore, they 
measured an adsorption energy of 1.0 eV. This value was 
obtained through an examination of the coverage depen- 
dence of thermally programmed desorption (TPD) ex- 
periments. Their work was also compared with several 
first principles approaches, including PBE-GGA^ but 
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FIG. 4. (Color online). Computed binding energy curves for 
azobenzene on Ag(lll). Experimental and theoretical results 
from Ref. 16 are shown. vdW-DF (this work) produces a 
BE in excellent agreement with experiment. The adsorption 
height is however overestimated by 16%, as discussed in the 
text. 
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dN-Ag to the experimental value and relax the system, 
the BE is reduced to 0.9 eV, and the aromatic rings tilt 
an angle of 6° relative to the surface. Thus, it appears 
that the interaction between the aromatic group is too 
repulsive within vdW-DF. 

One possible source of the height overestimate may 
stem from the fact that within the vdW-DF func- 
tional, the equilibrium lattice constant for bulk Ag is 
4.26 A (vdW-DF2 gives 4.32 A). This is larger than what 
PBE predicts (4.16 A), which is itself already too large 
compared to experiment (4.08 A). To test this hypothesis, 
we compute the vdW-DF equilibrium height and BE of 
azobenzene on a silver slab constructed with the smaller 
PBE lattice constant; we find that the adsorption height 
is 3.64 A, and the BE is 0.98 eV. Using a variant of vdW- 
DF which performs better at predicting the bulk lattice 
constant of Ag (optB86b+vdW-DF^, 4.11 A) results in 
poorer performance. Our computed BE and height with 
optB86b+vdW-DF are 1.54 eV and 2.85 A respectively. 
Similar results are obtained when applying the C09x ex- 
change functional with vdW-DF (C09x+ vdW-DF^ for 
BDA on the Au (111) surface. Here, the C09x+vdW- 
DF gives an excellent Au lattice constant (4.09 A), but 
overbinds by the same order of magnitude as obtB86b- 
vdW-DF. This highlights the delicate balance between 
short and long range forces which must be achieved when 
developing transferable functionals capable of accurately 
describing heterogeneous interfaces. 

In conclusion, we have used DFT calculations to inves- 
tigate the influence of London dispersion interactions on 
the adsorption of 1,4-benzenediamine (BDA) on Au(lll) 
and azobenzene on Ag(lll). A non-local vdW density 



functional (vdW-DF) was used to determine adsorption 
energies and the corresponding molecular configurations. 
In the BDA/Au(lll) system, we find that inclusion of 
vdW interactions produces adsorption energies consis- 
tent with experimental results, significantly enhancing 
molecule-substrate binding over PBE and stabilizing flat 
adsorption geometries. In the azobenzene/ Ag(l 11) sys- 
tem, vdW-DF results in better binding energies com- 
pared to other dispersion corrected functionals, albeit 
with some over estimation of adsorption height («16%). 
This work provides evidence for the relevance of the vdW- 
DF approach for the structure and stability of metal- 
molecule binding. 
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